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FOREWORD 


This report presents the results of work performed by- 
Lockheed's Huntsville Research & Engineering Center, under 
Contract NAS8 -21 146 for the Aero-Astrodynamics Laboratory 
of the NASA /Marshall Space Flight Center. * 


It is the last of a series of reports prepared by Lockheed 
from 1966 through 1969 under Contract NAS8-18036 and under the 
above mentioned contract. 

The NASA technical coordinator for the present study is 
Mr. Roger R. Burrows, S&E-AERO-GG. 
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SUMMARY 

This study considers the problem of flight scheduling in the planar 
two-burn minimum fuel rendezvous of an interceptor with a target vehicle. 

A set of equations is developed which takes position and velocity of the interceptor 
and of the target as input data. These equations allow calculation of five control 
parameters: the durations of the first burn, of the coast, and of the second burn, 
and the average thrust direction for each burn period. 

In order to set up the rendezvous conditions, Levi-Civita 1 s regularized 
variables and the corresponding orbital elements are used in contrast to most 
other papers on rendezvous problems. This brings several advantages compared 
to the use of polar coordinates; 

• An elliptic target orbit can readily be handled. 

• A near -circular coast orbit of the interceptor does not cause difficulties 
because the (badly defined) periapsis is not used. 

• The resulting equations are fairly simple and their numerical treatment 
is stable. 

The scheme uses some simplifying assumptions which, however, are 
satisfied in most practical cases: 

• The burn durations are assumed to be short compared to the duration 
of the coast. 

0 The Keplerian ellipses involved must have small eccentricities. 

® Instead of dealing with a variable thrust vector, the scheme uses a 
constant average value during each burn. 

Due to these simplifications, the scheme will furnish values of the control 
parameters which generally do not result in an exact rendezvous of the two vehicles. 
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However, updating these control parameters throughout the first burn yields 
a more accurate rendezvous after the coast and the second burn. 

Thu§, due to the simplicity of the equations (resulting in short computa- 
tion time) the scheme can also operate as a first-burn guidance scheme in the 
sense that the flight scheduling is done repeatedly based on current position 
and velocity data. 

In order to accomplish the rendezvous accurately by the second burn, 
a closed loop terminal guidance scheme based on measurements of the relative 
position and velocity is necessary. The approximations and simplifications in 
the present scheme are too rough for this purpose. 

Such a scheme, the Dual Phase Plane Method, has been derived and 
simulated in Ref. 3. Hence, in this report rendezvous missions are only simu- 
lated up to the end of the coast, where the terminal guidance scheme can take 
over. The second burn is handled without updating according to the latest 
values of the control parameters calculated at the end of the first burn. 

Simulation results for a typical rendezvous case are included in this 

report. 
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NOMENCLATURE 



D 


O' 


D- 


Aa., AB . 
J J 

6 



e 


E. 

J 

f 

I 


M 

O 



r 

s 


S 0’ S 2 
a 

T 


t 


semi -major axis 

regularized elements 

flight path angle 

total burn duration 

first and second burn duration 

element increments 

separation angle 

Krone eke r symbol 

exit velocity of the thruster 

perturbation function 

thrust acceleration 

interceptor 

earth's gravitational parameter 
earth's center, origin 
perturbing forces 
generalized forces 
distance from the origin 

regularized time, value for the interceptor at rendezvous 

increments of s during first, second burn 

value of s for the target at rendezvous 

target 

time 
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NOMENCLATURE (Continued ) 


coa 


tar 


u. 

3 


v 


w. 

J 


GO 


x. 

J 

X 

z 


coast duration 
mission duration 

"burn-out duration" = initial mass of I divided by mass flow rate 

I 

Levi-Civita's var iatle s 
velocity 

modified derivatives of the u 
frequency 

Cartesian coordinates 

thrust angle (measured from the fixed Xj -direction) 
side condition 


Ay-, G Q , G 2 , M, (p 


auxiliary quantities 


Subscripts 


j = 1,2 

The subscripts with the following meanings are mostly omitted in the preceding 
list. In the text; however, they are used in addition to the subscripts appearing 
already here. ° 

0 initial values of the interceptor, first burn 

1 coast ■ 

2 rendezvous, second burn 

3 initial values of the target 
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Section 1 
INTRODUCTION 

Optimum rendezvous problems are boundary value problems for differential 
equations (DEQ) where a certain cost function must be minimized by an appropriate 
choice of control functions. By the modern methods of the calculus of variations 
(for instance Pontryagiri's principle. Ref. 1) it is possible to solve problems of 

this kind exactly, but only with a considerable computational effort. 

s 

The goal of this study is to simplify the rendezvous' problem in an appropriate 
way, ' such that the results can be obtained within seconds by an onboard computer, 
but without losing too much accuracy (5 km position'e.rror ). 

The way to obtain simplifications is to introduce restrictions which are 
satisfied in most practical cases: 

© The interceptor's trajectory is assumed to lie in a narrow circular 
ring. 

Then the reauired velocity increments are rather small and can be attained by 
® short burn durations. 

This allows linearization with respect to these burn durations. 

« The thrust forces are considered as perturbing forces acting on the 
interceptor, I, whose unperturbed orbit is a Kepler ellipse. Only first 
order perturbations are considered. 

• During each of the short burns the thrust is put constant in magnitude 
and direction. 


This reduces the problem of calculus of variations to an ordinary minimum problem. 
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Even if all these simplifications are made the system of equations to be 
solved is quite complicated. Much depends upon the choice of the coordinates 
for describing the trajectory in the powered flight phases and upon the orbital 
elements used for characterizing the coast periods. Three possibilities are 
considered: 

1. Polar Coordinates Associated with Classical Orbital Elements 

This choice is striking because of the simple geometric meaning of the 
polar coordinates and the classical elements. Unfortunately this method fails 
in many cases we are concerned with, since the classical elements are badly 
defined for near -circular orbits (the perigee for a circular orbit is undefined). 

This approach has been the subject of earlier publications (Refs. 2, 3). In these 
reports a very efficient terminal guidance technique, the Dual Phase Plane Method 
due to I. Kliger and W. Trautwein, has also been described (see in particular 
Refs. 4, 5). 

2. Eevi-Civita 1 s Regularized Coordinates and the Corresponding Elements 

Although the application of this set of parameters yields more complicated 
equations, it is advantageous due to the "linearizing" effect of Eevi-Civita 1 s 
transformation (see Section 3). Transition through a circular' orbit causes no 
difficulties in these parameters. Most of the present report is concerned with 
the derivation of the control laws in this case. 

3. True Anomaly as Independent Variable 

The Kepler motion, described by direction unit vector and reciprocal 
distance as functions of the true anomaly satisfies a system of linear DEQ with 
constant coefficients (Ref. 6). Thus, applying these parameters has the advantage 
of Levi-Civita' s variables, yet transformations of that complexity are not used. 

Only a short description of the parameters and the corresponding DEQ 
will be given in Section 6. 
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Section 2 

THE RENDEZVOUS 

Let' T be a passive target vehicle moving on an elliptic orbit about the 
earth's center O. In the same orbital plane a steerable interceptor I is assumed 
to coast on an elliptic parking orbit. A rectangular coordinate system x, , x 7 
centered at O is used in the common orbital plane. The rendezvous mission 
consists of transferring I to T such that they meet with equal velocities and the 
least possible amount of fuel is used. The interceptor's engine is supposed to 
be ignited for the first time at a given time t = 0. 

The system of the two vehicles is characterized by the quantity 

f 1 = earth's gravitational parameter 

and by two parameters associated with the interceptor's thruster: 

T = initial mass of I divided by the engine's mass flow rate 
e = exit velocity. 

The situation at time t = 0 is given by the initial target data 
r = distance OT c 

j j 

v^ = initial target velocity 
7 ^ - target flight path angle 

and the initial interceptor data 

Tq = distance OIq 
Vq = initial interceptor velocity 
= interceptor flight path angle 
& — initial separation angle according to Fig. 1. 
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The rendezvous between I and T will be attempted by two burn periods (one 
at the beginning and one at the end of the mission) and a coast phase in between. 
This is the simplest strategy closer to reality than the two-impulse rendezvous. 

No variation of the thrust force is allowed; the engine of I is assumed to 
be either on or off. But since the mass of I decreases during the burns, the 
thrust acceleration f increases during the burns with time t according to 

f(t) = T ; ^ (2.1) 

where t* is the current total burn time. 

If the first and the second burn durations are denoted by Dq and 
respectively, the accelerations during these burns will be approximated by 



respectively. If the total burn duration 

D = D q + D 2 (2.3) 

is minimized, fuel optimality of the rendezvous (in our approximation)' is guaranteed. 

Since the optimal trajectory will be sought only among the ones with constant 
thrust direction during each burn, only discrete angles must be introduced in 
order to characterize the thrust. We choose the angle X^ between the thrust 
direction and the fixed x^ -direction at t = 0 and the angle X.^ with the same meaning 
at the rendezvous. 

Thus, the quantities Dq, D X^, X^ are the control parameters to be 
calculated from the initial data, while quantities s, a to be introduced later for 
characterizing the location of the rendezvous are merely unknowns in the mathe- 
matical problem. 
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Section 3 

EE VI - Cl VITA 1 S REGULARIZATION 


Regularizing is removing singularities from differential equations and 
their solutions by introducing new variables by appropriate transformations. 
Methods for doing this depend strongly upon the nature of the singularities that 
are to be regularized. In the case of the two-body problem in celestial mechanics 
the corresponding DEQ 

j2 

d x. x. 

2 1 + ^ -3 = °> j = 1, 2 (3.1) 

dt r 

have a singularity at the origin r = 0 , but in the solution x.(t) this singularity 
becomes manifest only when the vehicle collides with the central body. 

Regularizations of (3.1) have been known for a long time. In 1906 T. Levi- 
Civita (Ref. 7) found his regularization of the planar two-body problem, but only 
recently in 1965 it has been extended to three dimensions by P. Kustaanheimo and 
E. Stiefel (Refs. 8 , 9). The importance of these transformations lies in the fact 
that they produce not only regular but also linear DEQ for the Keplerian motion. 

In the sequel we give a brief outline of Levi-Civita' s regularization as 
well as a collection of the formulas we need for the further development. For 
the derivations we refer to Ref. 9. 

3.1 THE KEPLER MOTION 

Levi-Civita' s regularization consists of introducing the generalized 
coordinates Up u^ in the physical Xp x^-plane according to the conformal 
transformation 
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X| — Uj “ U 2 y — ' 2u ^ u 


1 2 


(3.2) 


and of introducing the parameter (regularized time) 

t 

s* = 




dr 

r(r) 


(3.3) 


as independent variable instead' of the time t, where r is the distance of the point 
(Xj, x^) from the origin and satisfies the relations 


2 2 2 2 
X 1 + x 2 = U 1 +u 2 


(3.4) 


The velocity components Xp x^ are transformed according to 


du j 
"ds* 


du, 

i 

His 


= T < U 1 X 1 + U 2 x 2 > 


(-U- X + u. x_ ) 


(3.5) 


The application of the transformation (3.2), (3.3) to the DEQ (3.1) of the 
unperturbed Kepler motion yields the linear system 


,2 

d u. 


ds" 


2 

4 + co T u. = 0, j = 1, 2 

2 I J J 


(3.6) 


2 . 


where is an energy constant given for instance by initial values r^, v of 
distance from the origin and velocity: 

2 

“i 2 = I < ^ - TT > (3 - 7) 
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Equations (3.6) are solved by 

u. = a. cos co_ s* + B . sin co T s* 

J 3 I J I • 

j = 1, 2 . (3.8) 

du. 

= COj {-a. sin 6) s* + (3 . sin 6)_ s*') 

J J 

The integration constants a., p . are referred to as the regularized elements of 

J <3 

the considered Kepler orbit because they characterize the orbit completely, as 
the classical elements do. 

3.2 PERTURBATIONS 


Equations (3.8) give a starting point for handling- the perturbed Kepler 
motion given by the differential equations 

,2 

d x. x. 

— 2 ^~ + M = p- (3-9) 

dt r ^ 


where p^ are the accelerations due to the perturbing forces. The generalized 
forces q^ corresponding to the coordinates Up u,, are 


q 2 = 2(Uj Pl +u 2 p 2 ) 

q 2 = 2 (“ u 2 P 1 + U 1 p 2 ^ 


(3.10) 


The presence of perturbing forces causes the frequency 0 ) to be variable (the 
value 6)j given by (3.7) being merely an initial value), and it generates an in- 
homogeneous term on the right-hand side of (3.6). Furthermore the elements 
a., p . defined by (3.8) are now functions of s* rather than constants. 

J J 

It can be shown (Ref. 9) that by introducing an independent variable, s 
by a transformation slightly different from (3.3), one can come up with regularized 
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equations where the frequency a>j is again constant according to (3.7). The new 
time transformation involves the semi-major axis a of the ellipse osculating the 
perturbed Kepler orbit: 

1 


r fi 

where v is the interceptor’s current velocity. When primes denote differentiation 
with respect to s, this parameter is defined by the differential equation 


t' = 



(3.12) 


where 


a 


I ” 



is the initial value of the semi -major axis. 


(3.13) 


Applying the time transformation (3.12) together with the conformal 
mapping (3.2) onto the differential equation (3.9) of the perturbed Kepler motion 
yields 


u." + 0) T 2 u. = F. , 
J 1 J J 


with the perturbation functions 


j 


1,2 


F. 

J 



a 






Ik V ) 


(3.14) 


(3.15) 
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The velocity transformation (3.5) reads now 


L l' 2 ^ a- (U 1 X 1 + U 2 X 2 


U 2‘ ' 2 ( " U 2 X 1 * U l x 2> * 


(3.16) 


Equations (3.15) are solved by the method of "varying the constant, " which yields 
equations 


u. = a.(s) cos (0 r s f 6 .(s) sin 0) T s 
J J I I 

u j' = j^~ a j( s ) sin w i s + P j( s ) cos S J 


j = 1, 2 


similar to (3.8). With the abbreviation 


w. = — 1 
J 


j = 1, 2 


(3.17) 


they read in matrix notation 


U 1 u 2 


W 1 W 2 


f cos a^s sm {Oj.s 

-sin OXj.s cos OJ^s 


a 1 (s) « 2 (s) 

3 x (s) p 2 (s) 


(3.18) 


The elements a.(s), (3 .(s) are now obtained by integrating 
J 1 


a. 1 = - — F. sin 60, s, 8 .* = — F. cos 0) T s 

J Ur J I j 1 


(3.19) 


The semi -major axis a can also be written in terms of the elements: 


1 2 2 2 2 

a = j (<x/+cx/ + p/ + p/) . 


(3.20) 
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Section 4 

DERIVATION OF THE CONTROL 'LAWS 

The equations yielding the guidance laws will be established according 
to the following ideas : 

• Describe all motions involved in terms of Levi-Civita' s variables 
and regularized orbital elements. 

• When two trajectories match in position and velocity at a certain 
time, their elements agree. 

In the sequel the subscripts 0, 1, 2, 3 (second subscript if there are two) denote 
values associated with the parking orbit or first burn, coast, second burn or 
rendezvous, target orbit, respectively. 

4. 1 MOTION OF THE TARGET 


The semi-major axis a 2 and frequency^ associated with the target orbit 
are given by (3.11) and (3.7): 







(4.1) 


According to Fig. 1 the target's initial position is defined by the coordinates 


X 13 ~ r 3 COS 5 * x 23 = r 3 Sin 5 


Inverting Levi-Civita' s transformation (3.2) yields the target's first two regularized 
elements; 


a 


13 


= ^3 


cos 


j6 

2 


a 


23 


\p3 


sin 2 


(4.2) 
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On the other hand, the initial velocity of T is 


"13 


= v 3 sin (-y 3 - 6), x 23 = v 3 cos (y 3 - 5 ) 


(4.3) 


as seen in Fig. 1. From the definition (3.17) together with (3.5) we obtain the 
remaining orbital elements 


P 13 
13 23 


200 . 


2oo^ 


( a 13 x 13 + a 23 x 23 ^ 


< -a 23 X 13 + a 13 x 23 ^ 


or by using (4.2) and (43) , 


#3 


3 V 3 


13 


23 


2oo, 


V F 3 v 3 


s in (y 3 - ~2 ) 


(4.4) 


20^ COS 'T'S - 2 > 


The quantities a- 3 , (3 ^ 3 are the initial values for solving the unperturbed regularized 
system (3.6) defining the target's motion. 

We now assign the values s = 0 and s = a to the initial point T 3 and to 
the rendezvous point on the target orbit, respectively. 

According to (3.18) is then given by the regularized coordinates 


u 12 U 22 


W 12 W 22. 


cos oo^cr sm u 


-sin 00- a cos a 


a 13 a 23 


13 p 23 , 


(4.5) 
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The time t^ taken by the target to move from T 3 to T 2 is calculated 


from (3.12), putting a/a Q = 1 and using (3.4) and (3.18): 


or 


cr 

' / [‘ 


( a, , cos GO, s + p sin 00, s)* 


tar / L'“13 2 ' h 13 ~2 

0 


+ (“23 cos s + [3,^ sin oo, s) ds , 


23 2 


* _ „ „ , a l3 P 13 + a 23 13 23 „ 

‘tar " a 2 CT + 2^ (i “ COS 2 “2 a > 

2 , 2 2 „ 2 
a 13 + a 23 " P 13 ' P 23 . „ 

+ is; sin 2 ^2 a • 


(4-6) 


Here we have used the relation (3.20) which can also be written as 


1 , 2 , 2 2 , 2 
a 2 = 2 u 12 +u 22 +w 12 +w 22 > 


(4.7) 


according to (3.18). 

4.2 THE BURN PERIODS 

In order to deal with the interceptor's motion we must define the 
quantitie s 


2 



analogous to (4.1). 


(4.8) 
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The motion of 1 will be described in terms of the regularized elements 
a j’ ^ j' Their values a^ 0 , (3j Q at the beginning of the mission (t = s = 0) agree 
with the corresponding values u. Q , w Q of the regularized coordinates. These 
are found by a sequence of transformations similar to (4.2) through (4.4). Taking 
into account the particular situation of Fig. 1, we obtain 


U 10 a 10 ^ 


v. 


w 


10 


= P 




10 


CO, 


u 


20 


= a 


20 


= 0 


(4-9) 


v. 


sin yo , w 20 = P 20 = ~ 2 ~ 


0 

^ COS ^0 


The intermediate elements (the constant values of the elements during the coast) 

will be denoted by a.,, (3 and the final elements a. OJ 6 (at the rendezvous 

JJ- j 1 j2 j2 

t 2 ) are related to the regularized coordinates u^, w^ 2 by means of (3.18): 


u 12 

U 22 \ 


/ COS COq s 

sm 0 Oq 

w 12 

W 22 / 


\-sin G)q s 

COS CDq 


/ “12 a 22 \ 

( 1 . (4.10) 

V P 12 $ 22 / 


The unknown, quantity s now stands for the value of the regularized time on the 
interceptor's orbit at T^- 


We further introduce the element increments 


Aa. - a._ - O' 

J J2 JO 


A(3 . = p - p 

H j2 M j0 


j = 1, 2 


(4.11) 


which are caused by the perturbing effect of the thrust during the burn periods. 


The next step is to relate the element increments to the control parameters. 
These relations are based upon an approximate solution of the DEQ (3.19) for the 
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burns. The right-hand sides will be evaluated at the rendezvous during the 
entire second burn (at the beginning t ~ s = 0 during the first burn). This is 
the principle of first order perturbations. Thus we obtain 


“ji 


"jO =° 


p jl " P j0 


“j2 ' “jl 


_ -p 


— F._ sin s 
“ft J 2 0 


j = h 2 


(4.12) 


P j2 " P jl 


— F.~ cos s 
“ n J 2 0 


where s^ and s ^ axe the increments of the regularized time in the first and second 
burn, respectively, and F. n and F^ are the values of F. at the beginning and at 


J 


the rendezvous, respectively. 


3 


From (4.11) and (4.12) there follows 


Aa. 

J 



— F j2 sin 0 ) 0 . 

S 2 S 0 
— F.-. cos s + — F.» 
j2 0 (£> Q jO 


(0 


0 


(4.13) 


and, by eliminating F 


j2 


from these two equations we obtain 





sin OOq s 




where 


(4.14) 



Aa. 

J 


COS (Op s + 




sin s 


(4.15) 


Equations (4.13) and (4.14) will further be used; to this end we prepare the expres- 
sion for F^ from (3.16); 
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r 2 q j + 




( q i V + q 2 U 2' 



(r 5 jk + w j2 w k2> 


*k2 


(4.16)' 


Here 6^ is the Kronecker symbol, and r^, q^ are gl ve]a hy 


2 2 

r 2 " u 12 + U 22 ’ 


(4.17) 


/’ 12 I 

■ 

2e 


^ U 12 

u 22 ^ 


1 cos A" 2 1 

1 q 22 1 


T - D o- 

D 2 

2 

l _U 22 

U 12 j 


i siaX Z / 


(4.18) 


according to (3.4), (3.10), (2.2). 
Thus, (4.13) is of the form 


The quantities u 


j 2 ’ 


w._ are 
j2 


defined in (4.5). 


r— 1 

d 

< 


/ cos X 2 1 

\ Aa 2 

= - M 

\ sin X 2 ( 


(4.19) 


where M is the matrix 


e • sin s 

“o 

/ a. 2 

f r 2 +W 12 

W 12 W 22 ) 


I u 12 

u 22l 

°2 

T " D 0 " 2 

2 * 

“2 

1 W 12 w 22 

. 2, 
r 2 + w 22 1 


i" u 22 

U 12 1 


The two equations (4.19) easily allow the elimination of the unknown thrust direc 
tion X ^ (we intend to keep only the unknowns s and o in the equations): 


(4.20) 
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( Aa x 


Aa 2 ) (MM T ) 


-1 


A “i \ = i 

A “2/ 


(4.21) 


where M is the transpose of M. We record a few intermediate results of the 
evaluation of this matrix product: From (4.20) we obtain 


, S 2 “o 

MM = r 2 ( I" —T ' 

“2 T - D 


-dT Sin “o s ) 


0 2 


,r 2 +w 12 W 12 W 22 


w 12 W 22 r 2 + w 22 


by using (4.17), and the inversion yields 


T -I 
(MM ) 


( s 2 a 2 r 2 


3/2 “0 


CO- 


D 2 

T “ D o - T 


sin COq s) 


(4.22) 


x 2 

2 + W 22 

-W 12 w 22 

12 W 22 

r 12 +W 1 


when (4; 7) is applied. Multiplying the matrix in (4.22) from the left and from the 
right by the vector (Aa^, Aa^) further yields the expression 

2 2 2 2 
r 2 ( Aa l + Aa 2 ) + ( r 2 + a 2^ ( w 22 Aa l " w 12 Aa 2 ^ » (4.23) 

when again (4.7) is used. 
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Finally, the time equation (3.12) is integrated approximately (first order 
perturbations) resulting in 


D o “ r 0 S 0 ’ 


D, 


00, 


00 , 


r 2 s 2 


(4.24) 


for the burn, durations Dq, D^. Using equations (4.22) through (4.24) in (4.21) 
now yields the equation 


d 2 

sin % S ' D 2 = G 2 " D 0 " ~Z~ ^ 


(4.25) 


with 


G 2 " “2 ^ 


r 2 (Aot 1 2 + Ao^ 2 ) + (1 + ) ( w 22 Aa l “ W 12 ( 4 - 2 ^) 


which contains no other unknowns than s, cr, Dq, D 2 • 

A similar equation can be derived from (4.14) in the same way; the 
result is 


*T sin “o s ■ D o 




(4.27) 


with 


G 0 = 'V ^ r 0 + A ^2 2 ) + <» + ^ > < w 20 A U - W 10 < 4 ' 28 > 

and A-y^ from (4.15). 

In order to obtain the total burn duration D we first solve (4.27) for Dq: 


D 


0 “ 


T G, 


A + 2 G 0 


(4.29) 
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where A is the abbreviation 


A = ^ ju e sin COq s 


(4.30) 


Using this in (4.25) then yields 


D 


^ T G 2 <A -j G q ) 

2 (a+Ig 0 )(A4-1g 2 ) 


(4.31) 


Hence the total burn duration D = Dq + is the function 


D(s, a ) = 


x(g 0 + g 2 ) 


(X Gq) (A+“ G 2 ) 


= min 


(4.32) 


which is to minimize according to the requirement of fuel optimality by appropriate 
choice of s and a . 


4.3 


THE COAST 


The coast trajectory of the interceptor is characterized by the intermediate 

fcs 

of (4.12) 


elements a^, 3 which can be obtained from the first and the last two equations 


a., = a. n 
jl JO 

j = 1, 2 (4.33) 

3 - n = 3 + Act. • cotan on s 

These quantities allow us to define the coast semi -major axis a^ and the cor- 
responding frequency 00^: 
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2 ( a l 1 + a 21 + P 11 + P 21 ) 


03, 




(4.34) 


JL 


4a. 


Now we can calculate the time t the interceptor takes for the coast 

coa 

by integrating {3.12) from to s - s^: 

s -s. 


coa 


^ja Q J ^(a 21 cos COj s + (3 12 sin C0j s)' 


+ (a 21 cos 0) JL s + p 21 sin cOj s) 2 J ds 


or 


coa 


= ^~ h (S ' S 0 ' S 2> 

+ a nh l + °h f> 2 l 


20 ), 


(cos 20 ^ s Q - cos 2o) x (s - s ) ) 


(4.35) 


a ll 2 +a 21 2 - Pn 2 - ft 21 2 
4a), 


(sin 2£0 1 (s - s 2 ) - sin 2a> 1 s Q )j, 


expressions for s Q , s ? _ in terms of s, a are obtained from ( 4.24 h 




^2 ^2 
(0 n r„ 


(4.36) 


The condition 


Z (s, cr) = t + D - t = 0 
coa tar 


(4.37) 


consequently guarantees that the two vehicle 
simultaneously. 


s arrive at the rendezvous location 
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Thus the rendezvous problem is reduced to the problem of minimizing 
the function D(s, a ) while the side condition Z(s, o ) - 0 must be satisfied. 


If once s and a are calculated, the burn durations are found from (4.29) 
and (4.31). Starting with (4.19) we will finally establish equations for the thrust 
angles X ' Q , X Equation (4.19) can be written as 


I Aon 


l Act. 


= const 


' r 2 +W 12 W 12 W 22 ' 


' W 12 W 22 r 2 + W 22 1 


l cos (x 0 - <P?) 


sin (X^ - <P Z ) 


(4.38) 


where 


<P 2 = aig(u 12 + iu 22 ). 


(4.39) 


Inverting (4.38) and forming the quotient yields 


tan {X 2 - <P 2 ) 


r 2 Aa 2 " W 12 * w 22 A “l " w 12 ^ a 2 ^ 
r 2 Actj + w 22 (w 22 Aa 1 - w 12 A« 2 ) 


(4.40) 


Starting from (4.14) we similarly obtain (since <P ^ = 0) 

ta„ V ~ r ° AT 2 ' W 10 (W Z0 At 1 ' W 10 Ay 2 1 
0" r 0 Ay l +w 20 (w 20 ^1 - w 10 A '<'Z ) 


(4.41) 
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Section 5 

THE CONTROL LAWS 


Here we summarize the equations derived in the last section in an order 
appropriate for computer programming. The numbers at the right-hand side 
refer to the corresponding equations in the previous sections. 


Input Variables 


/l, r, e, 6, r Q , v Q , y Q , r 3 , v 3 , y 3 


Constants (Independent of the main unknowns s, q } 


a o = 


0 


u 

2 


, v -1 

*2 * <4 ' 1 


CO 


C0- 


o = { ^ 


= yu/4z 


(4.8) 


(4.1) 


a 


10 = * 


a 20 = ° 


(4.9) 


V r 0 V 0 . y 

10 “ (O 0 2 Sin y 0 ’ 


0 


20 CO- 


U a, 

X cos y o 


(4.9) 


a 


13 


- fT 


cos 


a 


23 


i — . o 

- •yr-. sm 


(4.2) 




13 


co- 


V 3 y d . _ ^3 

sin ( \ ** ? )> P 

3 2 


v 


cos ( r 3 - | ) (4.4) 
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The Functions D(s,g), Z(s ,a) 


r u 12 

U 22 

] - 

1 COS M 2 <7 

sin a> 2 <7 \ 


j a 13 

a 23 \ 

w 12 

w 22 

I ' 

^-sin CC^CT 

cos 60 2 <7 


^ ^ 13 

^ 23 

r a 12 

a 22 


^ cos 0)q s 

-sin s i 

j 

! u 12 

U 22 

^ 12 

P 22 


' sin 00^ s 

cos 0)q s 

I 

1 w 12 

W 22 


r 2 U 12 +u 22 


Aa. = a._ - a.„ 
3 J 2 J° 


j = 1, 2 


Ay. = /la. COS a> 0 s + {(3 j2 - p jQ ) sin a) Q 


G 0 = “ 


2 , a "if 2 


0 ^ r 0 ^^1 + 4> 2 ^ + ^ + r 0 ^ ^20 A \ ~ ^ 10 ^ 2 * 


G 2 = “2 ^ r 2 ^ a l + Aa Z ) + {!+— ) (w 22 4a 1 - w 12 4ct 2 )' 


‘ X = | e sin fflp s 


(4.5) 

(4.10) 

(4.17) 

(4.11) 

(4.15) 

(4.28) 

(4.26) 

(4.30) 


D 


0 


tG, 


X + 2 G 0 


(4.29) 
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rA<G + G 2 ) 

D{s, a ) = D = — ~ x i (4.32) 

(A+^Gq) (A+l c 2 ) 



D 2 

= D - 

D o 


(2.3) 

Pjl ■ 

3 - t + ^ct. cotan 
H J 2 3 

a> 0 s 


j = 1, 2 

(4.33) 

a l 

2 (r 0 + P 1 1 + P 21 ^ 

®1 = 

: \h/Aa x 

(4.34) 

s o = ■ 

D 0 

r ’ S 2 

r o z 

“2 

®o 

D 2 

r 2 


(4.36) 

00„ 






o 

0)j ■ 

a l (s - s Q - s 2 ) 





+ !. 
20i 1 

( a 10 P li + a 20 P 

2 ^) (cos 2co^ 

Sq - cos 2oo^ (s - s 2 ) ) 

(4.35) 

+ 1 
2co 1 

( r 0 “ a *) ( sin 2co i 

(s-s 2 ) 

- sin 

2(£>i s 0 )] 


= a 2 ° 

+ 2&) 2 (a l3 P 13 + a 23 P 23 

) (1 - cos 2 {jl> 2 a ) 







(4.6) 

+ -^ r ~ 
2(0 2 

(r^ - a 2 ) sin 2C0 

2 a 





Z (s, (t) = Z - 

= t 

coa 

+ D 

*tar 

(4.37) 


The minimum problem D - minimum with the side condition Z = 0 can 
be solved iteratively without using derivatives of the function D(s, a ) for instance 
by the method described at the end of this section. 
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Output 


The variables s, o , s Q , s 2 , D Q , D D, t t are available from the 
computation of the functions D(s,cr), Z(s, a). In addition X Q and are 
obtained from 


tan 


r Q Ar 2 PjQ ^20 A7 1 “ P 10 
r 0 ^1 + ^20 ^ 20 ^1 " ^ 10 ^ 1 ^ 


(4.41) 


^2 ar § ^ U 12 + 1 u 22 ^ 


(4.39) 


tan (X 2 - <P 2 ) 


r 2 ^ g 2 ~ w 12 (w 22 Aa l ~ w 12 Aa Z ] 
r 2 4a 1 +w 22 (w 22 - w ]2 4« 2 ) 


(4.40) 


The Minimization Technique 

The idea is to seek the smallest value of the function D along the line 
Z(s,cr) = 0 of the (s, cr)-plane. A two-stage iteration process is used to improve 
an appropriate initial guess Sq, Oq almost to machine accuracy. A step size h 
indicating the order of magnitude of the error in the initial guess must be known. 


In the first stage for the three fixed abscissas Sq, = s^ + h, s 2 = Sq - h 

the one -dimensional secant method (starting from CTq) is used to find values 

= cr.* (j=0, 1, 2) which approximately zero the function Z(s, ct) at s = s.: 

3 1 


max 

j 


Z(s a.*) 

J J 


z < s o- V • 


(5.1) 


where is a small positive number, for instance - 0.005. 

A necessary condition for the feasibility of these operations is 


dZ 

du 


0 


(5.2) 
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in the region in which the arguments vary during the process. If (5.2) is violated 
the method works when the roles of s and a are exchanged. 


In the second stage the values D. = D(s., a.*) are calculated. Quadratic 

3 3 3 

interpolation then yields the approximate abscissa s of an extreme value of 
D(s, cr) with Z(s, gt) - 0; 


s 

m 


h 

2 


D 1 - D 2 
D 1 " 2D 0 + D 2 


(5.3) 


Finally, the corresponding value is calculated by quadratic interpolation 
with the collocation points and the values : 


cr.* - cr-* - o.* - 2o n * + cr-* 

a = or-* + (s - S n ) — — + (s - s -) — ^ ~ (5.4) 

m 0 ' m 0 2h ™ ?. < ' ' 


m 0 J 


2h 


The iteration cycle is closed by assigning the values s^, to the 


variables s^, Oq and by taking 


h = £ 2 < s m - s 0> 

as the new step size, where e 2 is another small number, for example 

€2 - 0.1. Figure 2 illustrates the meaning of the various quantities introduced 

here. 


Applied to the case of the functions D(s,o), Z(s,o) this technique con- 
verges very fast. 

Tabulation of D and Z for typical rendezvous situations has shown that 
the solution of the minimum problem is unique in the rectangle 

0<W Q s<^ , 0<co 2 o<^ 
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Fig. 2 - Minimization Technique 

Furthermore the line Z = 0 lies inside a narrow strip around the line 
<x>qS = and the partial derivatives of Z in this strip are quite large. Thus, 

the line Z = 0 is well defined. 

With the initial guesses 

s - Si- a - — ZL. 

0 2 a ) 0 0 “ 20) 2 

convergence was achieved in all cases considered. Results of a five digits' 
accuracy were obtained with less than 25 evaluations of the functions D and Z. 
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A Fortran IV subprogram FCT(S, SIG) for the evaluation of the functions 
D(s,o) and Z(s, a) consists of about 40 statements. A calling program (the 
program CONTRL) which determines the control parameters from the current 
position and velocity data, can be written with some 80 statements. The run 
time of the program CONTRL on an IBM7094 computer is in the order of 
0.1 sec. 
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Section 6 

OTHER PARAMETER SYSTEMS 


The set of the control equations recorded in Section 5 is not too complicated. 
However, many simplifying assumptions have been used. The most serious one 
is the linear approximation during the burn periods. The solution = y(h) of 
the system of DEQ 


'= g (7, s) , y(0) = y Q (6.1) 

was approximated by 

Yi = y 0 + 0) • (6.2) 


A much better approximation would be the value obtained from applying the 
trapezoidal rule to the DEQ (6.1) 


*1 = ^0 


4 * 
2 


g (y 0 » 0) + g (y^ , h) 


1 * 


(6.3) 


Here one would introduce two more unknown control parameteters, namely the 
thrust directions at the beginning and at the end of the coast. This would weaken 
the assumptions of constant thrust directions, but one would have to solve a 
minimum problem in four variables. 


This idea might be applied in connection with the third parameter system 

mentioned in Section 1. There is a good chance that the simplicity of the equa- 

* 

tions corresponding to these parameters compensates somewhat for the complica- 
tion of introducing two more unknowns. A summary of the most important rela- 
tions and properties associated with these parameters is given here. For more 
details see Ref. 6. 
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We start with the DEQ (3.9) of the perturbed Kepler motion and introduce 
the direction unit vector 


Yi = 


x. 

» j = 1» 2 


(6.4) 


and the reciprocal distance 


P = 


(6.5) 


These quantities define the position of the vehicle uniquely. If the parameter CD is 
introduced by 


dt = r df, t = I r d<p , 


= f 1 


d cp 


( 6 . 6 ) 


the dependent variables y^, p and t satisfy the DEQ 


5" + " iy j = ^ [ p j -5 E p k y k 


P" + Pf P = P - -T 


^ E 




, j = 1, 2 


(6.7) 


t' = 


(6.8) 


where & is the semi-latus rectum of the osculating Kepler orbit. S. is defined by 


1 ~ n ' ( X 1 dt 


dx 2 dx 

- X„ -TT- ) 


2 dt 


and satisfies the DEQ 


V = 


P P‘ 


• E p fc y k ' 


(6.9) 
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If there are no perturbations (p^. •= 0) Eqs. (6.7) are linear DEQ in y., p 
with constant coefficients (because of (6.9)). They describe a harmonic oscil- 
lator with the center y. = 0, p = y . Hence, a simple treatment of first order 
perturbations is possible, which is analogous to the method applied to Levi- 

Civita's variables in Section 4. In addition the parameters y., p and the cor- 
» \ 3 

responding orbital elements (Ref. 6) do not show any singular behavior in a . 

transition through a circular orbit. 

- In the derivation of the rendezvous conditions one can take advantage of 
the simple geometric relations given by (6.4) (6.5) and of the fact that in the un- 
perturbed case the independent variable (p is proportional to the vehicle's true 
anomaly. 
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Section 7 

COMPUTER SIMULATION 

Simulating a rendezvous consists of imitating on a computer all opera- 
tions influencing the trajectories of the two vehicles. The motion on their 
orbits is represented by theoretical or numerical solutions of the corresponding 
differential equations. 

For handling the coast phases it is necessary to calculate the position 
and velocity vector of the interceptor where it is Influenced by the earth's 
gravitation only. This is the initial value problem of the Keplerian motion. 
Levi-Civita' s variables allow solution to this problem in a stable and efficient 
way. 


From the vehicle's initial coordinates Xp x~, and initial velocity com- 
ponents Xp x^ (at time t=0, relative to an inertial coordinate system centered 
at the earth's center) the corresponding regularized coordinates and elements 
can be calculated by formulas of Sections 3 and 4: 


r 




a = 



(3.11) 


00 



(//is the earth's gravitational parameter) 


(3.7) 
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a 2 ~ 2a , 


^1 "" 2 to (Q 'l x i :f a 2 x 2^ &2 ~ 2oo (“ a 2 x l + a l x 2^ 

These quantities being known, it is possible to establish the equation 


At = acr + 


sinoocr 


(r-a) cosax7 + *i + a 2 sinoocr j 


which relates the true increment At with the parameter value o corresponding 
to the vehicle's position at time t = At. This is essentially Kepler's equation; 
it is most efficiently solved for cr by Newton -Raphs on' s iteration starting with 
the initial approximation 


"0 = a At 


Ml + M2 

2(0 


Finally, the transformations 


Uj U 2 \ 


1 w 2 / 


coscoa smcocr 


-smcocr coscoa 


a , 

1 2 


8 B 
Cl 2 


X 1 = u l 2 ‘ u 2 2 ’ x 2 = 2u l 
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• 2(0 , , 

X 1 = 2 2 (U 1 W 1 " U 2 W 2> 

U 1 +u 2 


* 2(0 , , , 
x 2 = (u 2 Wj +UjW 2 ) 

U 1 + u 2 

• • 

yield quantities Xp x^, Xp x^ representing now the position and velocity of 
the vehicle at time t = At. A subprogram KEPLER collects the set of the 
above formulas . 

The trajectory of the interceptor during the burn periods is calculated 
by numerical integration (for example Runge -Kutta) of the differential equa- 
tions 


x. 


V 


3 X 1 + 


T-t* 


cos X 


X 2 = 


1 1 , e 

3 x 2 + T-t* 


sin X 


(7.1) 


in the Cartesian coordinates Xp Xp where X is the current thrust direction 
measured from the x-^-axis. 


The two computer programs CONTRL and KEPLER as well as the sub- 
program INTGRT for the Runge -Kutta integration of (7.1) are put together to 
a simulation program according to the rough flow chart shown in Fig. 3. 


The control parameters are updated during the first burn only, when the 
logical variable UPD is TRUE. The fixed updating time interval is DT, while 
TINT denotes the current updating interval. 


The computer output sheets (Fig. 3) give an account of the two vehicles' 
motions and of the control parameters for two test cases. The lines in the 
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INPUT 


Earth's gravitational 
parameter 

V 


Burn-out 

duration 

r 


exit 

velocity 

e 


Updating 

interval 

DT 


radius 


velocity 


■ flight path angle polar angle 




>3 

^0 


6 

<P 



Fig. 3 - Flow Chart of the Simulation Program 


target 

interceptor 
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Fig. 3 


Flow Chart of the Simulation Program (Continued) 










LMSC/HREC D149128 


section Control correspond to successive iterations in the minimization procedure. 
The last line gives the values of the control parameters. 

As an example, the circular rendezvous case mentioned in Ref. 2, page 9 
is considered. Initially, the interceptor and target are on circular orbits 100 km 
and 400 km above the earth's surface, the target being 5° ahead of the interceptor. 
The table in Fig. 4 compares the results of the present simulation program with 
the exact calculus -of -variations solution (COV) and with the results of the guidance 
scheme in Ref. 2. In the present scheme three cases are considered: (1) no up- 
dating (DT =«>), (2) updating interval DT =10 sec, (3) DT = 4 sec. 



COV 

Present Scheme 


Ref. 2 



DT= oo 

DT=10 

DT=4 


First burn Dq (sec) 

13. 29 

13.20 

13.29 

13. 31 

13. 3 

Coast t (sec) 

coa 

2522. 7 

2505. 2 

2512. 2 

2508. 9 

2318.4 

Second burn (sec) 

11.97 

11. 76 

11.96 

12. 00 

12. 1 

Thrust direction X q (deg) 

64. 9° 

64. 6° 

64. 6° 

64. 6° 

62.8° 

Thrust direction at end of 
First Burn X ^ (deg) 

66. 1° 

64. 6° 

66.3° 

67. 0° 


Final thrust direction X^ (deg) 

-129. 4° 

-106. 0° 

-106. 5° 

-105. 0° 

o 

CO 

CO 

rH 

r—H 

1 

Position error (m) 

0 

5616 

971 

566 


Velocity error (m/sec) 

0 

98 

100 

100 



Fig. 4 - Comparison of Simulation Results 

In the second example, the results of a rendezvous with a target moving 
on an elliptic orbit with semi-major axis 6793 km and eccentricity 0.0100 are 
shown. The interceptor is initially 8.6° behind the target and 6400 km away 
from the earth's center. Its orbit has the eccentricity 0.0246. Initially, the 
first burn, the coast and the second burn are predicted to last for 35.0 sec, 
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2600.7 sec and 15.2 sec, respectively. Updating throughout the first burn 
in intervals DT = 10 sec modified these numbers to 35.8 sec, 2628.3 sec, 

15.9 sec. The final position and velocity errors of 9317 m and 137 m/sec, 
respectively, are in the same order of magnitude as the errors in the cir- 
cular case. The results in this case are collected in Fig- 5. 

Further testing showed that the present flight scheduling and guidance 
scheme yields good results when the assumptions of the scheme — small 
eccentricity in all Keplerian ellipses and short burn durations — are 
satisfied. A closed-loop guidance scheme being applied to the second burn, 
however, would necessarily use a burn duration somewhat longer than the pre- 
dicted one in order to compensate for the more pronounced velocity errors. 



COV 

Present Scheme j 



DT= oo 

DT=20 

o 

i — t 

ii 

H 

Q 

First burn D^(sec) 


34. 96 

35. 65 

35. 82 

Coast t (sec) 

coa 


2600. 7 

2626. 2 

2628,3 

Second burn (sec) 


15. 16 

15. 64 

15. 90 

Thrust direction X^ (deg) 
Thrust direction at end of 


129, 8 ° 

129. 8° 

129. 8° 

First Burn X^ (deg) 


129. 8° 

132, 0° 

135. 8° 

Final thrust direction X ^ (deg) 


-83. 5° 

-81. 6° 

-80. 4° 

Position error (m) 

0 

37653 

29842 

9317 

Velocity error (m/sec) 

0 

144 

145 

137 


Fig, 5 - Elliptic Case 
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—0.1a 45 7 c 04 

-U. 753126 


HfcNutiVl Ub 











, 




bfci 'Mi* mUh- 


. 1 6 ] bPo I n- — ■ . 
B . 251tt a 5 - * 

2 . 25 J. 931- .1 

0.2520/b »•-* 
O. 2521 9b . a 
u. 25231b -4 
!■). 25242 b "a 
0.25254b ► a 
0 . 25263b " *i 
0,252/86 a 
B. 25269S r4 
0.253016 -a 


•■■AnlhSiAN i.Qi.R01~M»fcb 


ANl) 


JI . 6 6 o 2 6 c 07 
•1 • 60 32* 4 c _ 0 7 
2.60361b 07 
2.69-i,7 9£„l'7. 
4.63396c 0/ 

/ .66 A L4c 0" 
3 . 6 6 A 3 1 fc If. 7 
2 ■ 6p448 b . jo 7 
/'.6o465b 0 7 
• 6j>.Ad2c 0 7 
/ .63 A 76b 3/ 


0 . 13590b 0 7 
. 0., 1 35 03b 0 / 
l). I34l5t B7 
0. 13328b 07 
0. 13241b 07 
K.13l5ob 0/ 
0.13066b 07 
0.1 2978b 0 7 ■ 
0.12891b b7 
0^1 06 04b 0 7 
■0.12 716b' 0 7 


_yt. | .t , U.l Y 

0.l5l/2b 
0 . 1 5069b 
1 490&t 
148.63b 
14760b 
14656b 
14553b 
14450c 
i’43A6b 
1 4 ’2 4 3 b 


-0.14139b 


CO MPONENTS 
04 -0. /4292b 04 
04 -0.74302b 04 
>4 -0./43126 04 
04 -0. 7 4.3215. 04 
04 -0.74331c 04 
04 -0.743405 04 
04 -0.74349b 04 
04 -0 . 743585 04 
04 — 0. 7436/t M4 
04 -0. 74376b 04 
*•4 -0 . 74365b 0.4 


HtNub" ‘ 


(. 1 kS I . Sbt 1 Mi C»7N. MiSoI'-JN 4.l3l950b 02 fc.il/578b 02 

~ Rbi Al l Vt ► l 9 I t 1 lift VtC * .-4 4. 1951 03b «4 -0.526o50b «4 

Hbl.Allvb ■ri.IH.n* >7 2.317^4 04b 0? 0.9"74flb V'.2 

39 


0_. 253016*- _0_4 
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CIRCULAR CASE, UPDATING AFTER 10 SECONDS > 


Input 

MU ] ft u 

b 

m 

yi . 3V603bb 

lb 

0. 583326b 

03 

O.418/40E 04 

0.100000b 

H d y A 

_l-Aq_ 

UEL 

M.&/7530E 

0/ . 

0 . 766874b' 

04 

0. 000000E- 00 

0.872665b- -1 

P* 3 '/ W 

LAP 

PHi 

. 647o«$it)fc 

0/ 

0. 784442b 

04 

0.000000b 00 

" 0 .000000b /-0 


StM i - ',A JCiK AXfco_uf I N > b R CbP rOK , ..TARGET 0,64733b V>7 0 . 67/32E 07 ‘ 


lAPGfcJ t I Ml- 
b v fcl 

LAf<JfcSlA/>j UOORuI'^AT^b A NO 

tf. 6/4 75 t 7 tf.bVtf'S’Sb 156 

ybLOCirV COMPONENTS 
- 0 -. 668 o 7 b E 3 0 ./ 6393 E 04 



- 


LjN > Rut 
b 

b 1 u 

l D ■ 

l ' 

. 00 _LOa * 02 

CH 0 

CP 2 

. 40 # 48 b-fc«O 

0 . 4 i» 9 o 6 c-li 3 - 16 / 
3. 60 /OJ.C -03 -0 

. o 2 v'l 7 9.3602 

.0623 24 . 9720 ' 

42.3021 232 / 
13.2336 2462 

.7815 37.0580 

. 0/84 • 11. 7363 

-«.^410 
1 . ' 175 

-a . >.797 
- 1,0166 

4 . 37 /llb -03 

0 . 5 / 533 C-M 3 M 

.V» 0 f, H 24.9531 

13.1932 2503 

.2168 11.7578 

• . 1294 

-l.o 513 


LA^TbSlAfi CGOROlNATbS ANO 

(E'.OuliY COMPGNENIS 




0 . 6 / 4 w 4 b 07 

0 ,obb 68 t 06 

- 0 . 73482 b 03 

0 . 76513 b 04 




CAofbSUN COORuiNA fES AND 

VELOCIlY COMPONENTS 



a. 00 a 00 b Lb 
W, 10000 t 01 

4 . 64 /o 3 b 07 
7 . 64/328 t )7 

iJ . */ 0000 e 00 
0 .- 78476 H 04 

0 .M 0000 E' 00 
- 0 . 64366 b 01 

0 ./ 8444 E 04 
0.785096 04 



0 . 2020 Mb > 1 
W. 30 fO 0 b l-l 

B. 64 / 32 t 07 
7 . 64732 b 07 

0 ,l 5 / 01 b 05 
0 . 23362 b 05 

-tf, 128 t> 7 fc 152 
- 0 , 1 V 293 fc 02 

0 . 78574 E 04 
0 '. 78638 E 04 ' - 



/' . 0 9 y 9 9 b il 
U. 49 y 9 Vb J 1 
0 . 6 9 y 99 b K 1 
a. 6 yyv 9 b fc 1 

t , 6 ^/^ 2 fc U 7 

0 , 31429 b 05 
0 . 39303 c 03 1 

.- 0 . 23/148 02 
- 0 . 32129 b 02 

0 . 78 /W 3 E 04 
0 . 7 b 768 fc? 04 



v. 64 751 c 07 
0 . 64731 b 07 

0 . 47183 b 05 
" 0 . 55069 b 05 

- 0.383398 02 
-H-, 44943 b 02 

0 . 78635 E 04 
0 . 78898 b 04 



tf./ 9 y 99 b 01 
^.OVvVVP * 1 

4 . 64730 b 07 ’ 
0 . 64 / 60 C 07 

~ 0 .- 62 u 62 c- 05 • 
0 . 70862 b' 05 

. - 50 . 51342 b- 02 - 
- 0 , 3 / 735 b 32 

0 *. 7 . 8963 E JJ 4 
0 . 79028 E 04 , 


‘ 

0 . 9 Jyy 9 b 01 

3 . 64729 b ' 07 

0 . i 87688 05 

- 0 . 64125 b 32 

0 ./ 9093 E 04 




CJftS 1 -* JL . w 

^ b 1 v> 

0 . 3 / / llb.-M.j_ M .3/ 3- j5b-03 
0.579/9C-05 U • 5/ ?14 c->'j3 


. 7 

-36.5366 
Id . 0000 


0 i 00 JOA 02 CO0 CH2 

2b. 4900 10.o402 2 466. J. 4 1 1_ 14 ; 8 498_r /• • '<$5_6._-0. *J>£/. 

15.2585 0.2941 25lT.lt> b5~ 11 . 9<T4 4 1..583 -1.<j409 


! 

[Akbh i I lr.fr . UA ofcSIAN UOQRu 1 WA I bS AN p V. EL PC If Y COMPP N§HT S_ 

i».l3«94h ici ” /,o.o/6c 07 /)',o9l82b 06 , -0.78328E 03 0.762&6E 04 


lM.eKi, T I ML . 
’t0.y9y99b 01 
M.IMJbVfc 02 
0 . 10o56t 02 
0 . 10y8ob .02, 
0 . 11 o 1 / b 02 
0.11s4/£ P2 
0 . 11 y /6c 0'2 
_W . 1 2005b 0 2 
0 . 12fc3bb L2. 
0.12VO4C K'2 
0.13294C K2 


..•A. . cS I A:-) COOftOlNAffeS ANU 'vELUCirY COMPONENTS 

'.64/29C 07 0 . / 6 /66b 05 ' -0,64123b 02 0.79093E 04 

f . 64 729t_ 07 a. 61374 b 05 _jJ 4.6 6 269b. j>2_._.0;_7.91 .1 4E _04 

2, 6 fi "/ 29b 07 0.63980b 05 ' " -0 . 680SI '02 0.79136E 04 

'28t_02 M .6o3Q/ b 05" -0 .70 620b 02 0 . 7.915 /£ 04 

2.6W28E j7~"’'d,b919ob 06 -0.7^/o4t 32 0.79179b 04 

0 . 64 / 28b 07 0. 91604b 05 ' -0. 74948b 02 0. 79201c 04 

0.64/ 28b 07 0 . y 4 4 1 3b 05 -0.7/Ulb 02 0 . 79222b 04 

■A . 6 <» / 2 7 e 07 0 . 9 7023c 05 -0. 79274c 32 0 . / 9244b 04 

0.04727b 07 0.99003b' 05 - 0 . 6 1 4 3 6 E 02" 0./9266E 04 

4. 64727c 07 3. 10224b 06 . -0'.8.>397t 02. . « _• 7.9287b. 04 

"7.64/26C 07 *J.104&3b .06 """-0 . 63758c' 02' ‘ a. /9309b 04 


bNU Of ' f )KS I -oMh'' 


. 40 
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i AkUb l i 1 i*>t . 


CAm r t- b> X A iv 

6 1) OK U I iv A T hb 

A i\l -J 

VEU-iLlIY U-WU'YfcNTS 


0 . <t 8 o / ' •* r. » <♦ . 


- 3 . o o o 2 1 1 

07 

0. 1.2222b 

0 7 

-0.13858b 

0 4 -0. 75428b 

04 

Rfc^'Ub >. vuii6 

















- 

St(,U-\'.l 1- U rx ■' 









1 N 1 tb'v i 1 Mb , 


l. A n i b S 1 A M 

C 0 0 R u i "i A ! h b 

A N U 

VbLOUllY CUHPONEMS 

- 

6 . y a 


-/l.6o461b 

0 7 

0 . 13122c 

0 7 

-0.146oiE 

04 - 0 , 7 4 5 7 8 fc 

0 4 

0 . 2 8 2 0 6 b i' A 


-0.6o4o8fe 

0 7 

0 . 1 30o3i: 

0 7 

-0.14826b 

04 -0.74584b 

04-‘ 

0 . 2 8 £ / b b *- 4 


— Si . 60486b 

07 

0 . j 2y< : 4b 

0 7 

-0 , 14421b 

04 -vj. 74395b 

04 

0 . 28</90b 1 A 


~'i . 6o6vj3b 

07 

0 . .1 26 t 'Ob 

07 

-0. 143x6b 

04 -0.74402b 

04 

ir>.25oa2fc v 4 


-Si . 66 6 20c 

a 7 

0 . i 2 / r»6t 

07 

-0.14211b 

04 -0./4411E 

04 

0 . 2 8 O 1 A fc- ■' 4 


- S' . 6 o 8 5 7 b 

07 

0 . 1267 7b 

07 

-0 .14106b 

04 -0 . 74420b 

04 

0 . 28 0-iOt !£ 4 


-■/ . 6 o 6 6 4 b 

07 

0. 126b 8b 

07 

-0. 1400lfc 

04 -0. 74428b 

04 

V) . 2 8 o So b ('■ 4 


-/ . 666 / lb 

0 7 

0 . 1 2 4 <- V fc 

07 

-0 . lob46b 

04 -0.7443/b 

04 

0 . 25580c '-'4 


- vj . 6 o ? 6 7 b 

0 7 

0 . 12 4': 0b 

0 7 

-0.15741b 

04 -0.74448b 

04 

0 . 2 8 0 6 2 b t 4 


- Si . c o 6 Si 4 1 

0 7 

M . j. 2 5 c lb 

0-7 

-0 . loobob 

04 -0.74480b 

04 

0 . 2 8 5 / A b v. 4 


- 0 . 6 0 5 2 0 c 

07 

0.1 220 2b 

07 

-0.15880t 

04 -0.74461b 

04 

RfcNUfc^VUUb 


















b 1 R S 1 . b t C t: iv U 

,i[jKi\i . M l b 6 t UM 

0 . J o2 l - 40 k 

0 2 0 

.1] 9 6 4 o h 

02 0.285741: 

: 04 

RfcLAl L v tr HOb i l 

1 u \ 

VC.UI ' i S' 


O , 6801 00 r 

32 k 

. 969280b 

05 


Rbl. A r t V'b Vt l.Ou 1 1 y 

f oil s^i'lbNT 

s 

0 .26/ v- 16b 

0 2 0 

.-9 66 V 72b 

02 



















t 













































































ELLIPTIC CASE, UPDATING EVERY 10 SECONDS * 


LMSC/HREC D149128 - 


INPUT 


MU 

• 1 AU 

t 

U1 

w . 398335b 

TIT 7 

W ; 583326b 


0.418/4OE 04 

0. 10W0W0b 


R3 

■J 3 

>;a i 

DEL. 

W. 680W00E 

0 / 


04 

-0.99999VE-02 

0.150000b 

■0 

RW 

v 0 

t- AW 

PHI 

0.6 40000E 

07 

8. 780000b 

W4 

0 . 999999E-02 

0. 00000BE 

40 


StMI-fiAJOR AXES Oh I NT ERCEP T OR , TARGET 0.62592S 07 U.67934E 07 


TARGEI TIME, 
0.00000E 00 

cartesian cooruinates and 
0.6/236E 07 0.101616 07 

VELOCITY COMPONENTS 
-0.12187E 04-. 0,755226 04 





CONTRUL 

13 

Sic 

Z 0 

00 • COA ' 02 

CH0 

CH2 

0.39380£-»3 
0 .39522E-P3 

0.41027b-03 - 
3 . 3V254b-B3 

•254.63W9 145.8677 

-0.02*4 50.12*4 

90.5155 2368,9336 55.3522 

34.9603 2599,0747 15.1641 

2.8885 

2.2648 

-0.0603 

-1.4558 

B .40289E-03 
0.J9 532E-03 

J . 4 

0 . 3V265>t- J3 

0.1196 50.1757 

-0.0415 50.1275 

35.0400 2651,9336 15.1356 

34.9642 2599.7549 15.1633 

7.2687 

2.2650 

-1.5063 

-1.4562 

0 .39543E-03 
0.39546E-03 

rt . 3y2d0t£-tf 3 
w . 3y? / dh-^3 

-0.0010 50.1203 

-0.0005 50.1204 

3*. 9592 2600,8369 15.1611 

34.9590 2600.7295 15.1614 

2.2648 

2.2648 

-1 . *579 
-1.4578 


TARGEI Tint, 
0.100008 02 

CARTESIAN COORUlNATbS AND 

0.6/ll0fc 07 0.1H916E 07 

VELOCITY components 
-0.13038E 04 0./5389E 04 




1 N TbRC TlMt, 
0.00000R 00 

lARTESIAN COOR'J INA1 £s ANU 

V,6*000£ 0/ 0.W0B00E 00 

VELOCITY COMPONENTS 
0.7/998E 02 0.77V96E 04 



0.10000b 01 
0.20000b 01 

# . 6 * Vi 0 i £ 

0? 0. 760236 04 

07 0.15610E 05 

0.63678E 02 0.78051E 04 

0.49350E 02 0.78106E 04 



0.30W00E 01 
0.09999b 01 

Y . 64fok3lfc 
A . 

•07 0.23423b 05 

07 W.31242H 05 

0.35014E 02 0.78161E 04 

0.206/IE 02 0.78216E 04 



0.49999E 01 
U.59999E 01 

7.640W2E 
v’, . 64001 E 

07 0.J9066E 0b 

07 0.46896E 0"5 

0.63203E 01 0.78271E 04 

-0.80385E 01 0.78326E 04 



U.09V49E 01 
0. 799 99E 01 

c. 64003 b 
4.64001b 

07 0.54732E 05 

07 0.O2573E 05 

-0.22405E 02 0.78381E 04 

-0.367/9E 32 0.78436E 04 



W.B9yV9b 01 
0/. 99999E 01 

£ . 64001E 
0 . 6*WW0b 

07 0. 70419b 05 

07 0.78271E 05 

-0.5H62E 02 0.78491E 04 

-0.65552E 02 W.78546E 04 




. . 

CUNT RUL 
S 

SIC 

2 D 

•00 COA 02 

CH0 

CH2 

W.39546E-1-3 
0 .39B83E-P3 

0.392/8E-03 
W . 39362b-w3 

-33.8467 49.7338 

-0.0015 40.7866 

31.6734 256/, 2373 18.0604 

25.3889 2615.7324 - 15.3977 

2. ’>155 
2.2748 

-0.6844 

-1.4419 

0.39B66E-03 

M.39345b-03 

0.0005 40.7864 

25.3868 2614,5205 15.3996 

2.2747 

-1 ,'4408 

TARGEl HUE, 

CARTESIAN COORUINATES AND 

VELOCITY COMPONENTS 



0.20000b 02 

0.669/5E 

07 0.11669b 07 

-0.13888b 04 0.75246E 04 



INTbRC TIME, 

CARTESiAN COORDINATES AND 

VELOCITY COMPONENTS 



0.99999b 01 
0.10999E 02 
w.livyye 02 
0.12999E 02 

2 . 6*000b 
3.63999b 
0.63998b 
0.63997b 

07 0, 7827 IE 05 
07 0.86128E 05 
07 0.93991b 05 
07 0.101856 06 

-0.655528 02 0.78546E 04 
-0. 80005b 02 0.78601E 04 
-0.94467E 02 0.78655E 04 
-0.10893E 03 0.78709E 04 



U.13999E 02 
0 .149996 02 

i. 639966 
0.639956 

07 0.10973E 86 

07 0.117616 06 

-0.12341E 03 0.787646 04 

-0.13790E 03 0.78818E 04 



0.15999E M2 
0.16999E 02 

2.63993b 

0.63992E 

07 0.125496 06 

07 0.13338E 06 

-0 . 15239t 03 0.78873E P4 

-0.16689E 03 0.78927E 04 



0.17999E 02 
0.18999E 02 

0.63990b 
0 .63<>88b 

07 0.14128E 06 

07 0.14918b 06 

-0.18140b 03 0.78981E 04 

-0.19592E 03 0.790356 P4 



0.19999b 02 

2 . 63986E 

07 0.15708E 06 

-0.21045E 03 0.79090E 04 




CONTROL 

b 

V.S9bbbb-Y* 

bib 

0.3V345&-03 

7 0 

-35 . 87 80 41.5457 

U0 COA 02 

22.9310 2577.8579 18.6146 

CH0 

2.6451 

CH2 

- io . 4 3 P 

0.40i//b-t3 
0. 40201E-00 

W . 3v5vwe-03 

-0.0093 31.3362 

-?* . 4316 37. 33d9 

15.6655 2626.6616 15.0407 

2ki.3197 761^.1743 17,0192 

7 . 2976 
2.5654 

-1. .248 
-1. -105 

0. *01858-03 
__0_._40l91E.t03 

0 . 393926-03 
_0 ,3.93x8c-43 . 

-0.021 5 31.3088 

-7-.0_0.24 3 1 ,_30 5 0. 

15.0685 2627.1577 15.6404 

" - 1 o 6 6.2 _ 26 2J j, 5_ 6 7 9 ..1.6^0 3.94. 

2.2979 
_ i. ,.29.77.. 

-1 . *250 
..-a.G-.25 7_ 


42 


NOT REPRODUCI 


ry 

lx 


i 
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TARGfcT 7 1 rib , 0 t ,;7t S 1 h \i “J 004 u 1 N A • c s AMU vfuOCliY COMPONENTS 

yi.Q-QOMOr- 02 0.6b832b O'/ 0,12421b 07 -0 . 14736b . 34 0 , 75g93 , E- . g 4 


INTbRb 'lint, 
0. 1999 9b 1 2 
0 . 20999 b t 2 
o.aivyvb . a 


0.22999b 1-2 

... 0. 2 3999b 02 

0.24999b ',12 

._J)...2.L9.99.f.>J(£ 


0 . 2699 9b 02 
a. 27999b >2. 


CAKlbSIAN COORDlNATth' 
0.6596 6 b.. $.7 _4 jl IjJ.i.O.&b 
<■.639846 07 0.1,6500b 

2..639d H=_a7 o 

0 . 6-59 / 9t 0? 0.1«084fc 

0.639/fi t: fl7__Q^_lM76b 
0.6->9/4t 07 0,19670b 

0. 669 / It 07 jj_. 20463 6 

0.63967b 07 0,212566 

0 . 6396 4E 07 0.220536 


AMO VELOCITY CuhPONbNTS 

0 6 -0.210 4 56 03 0.790906 0-9_ 

06 ” -0.225116 

06 
06 

SLLl. 


06 


06 

06 


0.28999b I 2 0 

0_^9J7.9.V.h_ fl2 


63961b 07 0.22848b 

63 967b 07 0 .236 44b 


03 

-0 ■ 239 7 9b 03 


-0,25447b 03 
-0 . 269 l6.tL.03_ 


0.79143E 04 
'0.79196E 04 


0 . 79249E 04 


-0.28386E 03 
-0.29856E 0 3. 


0.79355E 04 
O.79408E 04 


-0. 31328E 03 
-0.32801E 03 


0 . 79461b 04 
0.79514E 04 


06 -0.34274E 03 0.79567E 04 

06 -0. 35748b ...03 0_-._Z9.8ZaE _04.. 


CUNTROL 

_S SJLt L B 

0.40l9lfc-/i3 0.39398E-03 -37.5737 33.9510 

. ,_0. .40445 1-0.3.. _.a*3.9 3_6Ah-0 3 - 0 ..0 0 1£L . . 2 3- 7.150- 

0 . 40350b -03 0,392/0b-w3 0.0010 21.7160 

o .403486-03 a-,.39^& iih7:.0.3 0-._2 D .afl 21 .7 1 6ft . 


ne C.OA 02 CHfl CE2 — 

14.8405 2587.3291 19.1105 2.9042 -B.8075 


“ft 81 6 7 _ 26 28'! 4 4 78 ^ft8994 

_5jJ>lfc4_2.fe£A.jll. 7 ,4 1 ? . 8 9-9 2_ 


2.3707 -1.-026 
2.3707 -1 . «025 


TARGbl TIME, CANlsSiAN COURDlNATbS AMD VELOCITY COMPONENTS^ 

jO , 35j41Afc_il2 fl.jL.6A7 4j.t_ a 7._ .0^.12355 IX— HJ rDjJt '>223 b_ 0.4 — 0. 7 5 OD D b — 0.4. 


INTbRU I 1 Mb , CAKlhSiAN COORDINATES AND 

O . 29 9 99b .0.2 j3_. 64 V 5 7b_ J0? ._ 0j.21o4.4fc_ 0JL 

0.30581b 7.2 3.60955b' 07 0.241O7E 06 

_HU.311. 6_Ab Y 2 9. aJ.ni 07.. . 0 . 2i57ob._.B6 .. 

0 . 31 / 44b 02 3.63951b 07 0.25034b 06 

0 . 32o26b o2 0.63949b 07 0.25498b 06 . 

0.32908b 02 O.o3946fc 07 0.259blb 06 

0.33489b 02 0.6 3944b 07 0. 26425b 06 

0.340/lb 02 3.6c44lt 07 0.26890b 06 

0 . 3 4652b 772 t ,609 39b 0 7 _0_..2 7 354b 06 

0.35234b 02 .'..63936b 07 0.2/8186 06 

0.35alob 02 '0.6o934b 07 0.28283b 06 


VELOCITY COMPONENTS 
-0.35 748b 03 0 1 , .7_9 _ fr2flg_M_ 

-0 . 3 o629b 03 0.79649E 04 

-0.3 7 5llb 03. O jl- 7 967 7 E M. 

-0. 38393b 03 0.797O6E 04 

-0.392756 03 0.79734E 04 

-0.4O157E 03 0.797636 04 

-0 .41040b 03 0.79791E 04 

‘-0 . 4l923b“o3 0.79820E 04 

- 0 . 4 280ob 03 0.79848 E 04 

“-0.43690b 03 0.79877E 04 

-0.44574b 03 0.79905b 04 


END Oh \ I HS 1 bUfiN 


~T ARGfc I " T I Mb 1 “uAprbSlAiV oOORljlNAIbS ' ANO ‘vElOUTY COMPONENTS 

0.26800b_7 _4 2.Z. ■ 6 7 640b 07 -0.46 714 b 0 6 0,45317b 03 -B.76590E 04_ 

KtNUfcflVOUS . — 


SbCONu dUKN 


i N T b R 0 lint, 

0.26641 b 0 4 

0. 2665 7b 04 
0.26673E 04 
0 . 26eavh 04 

_0_. 26/046 _1_4 

0.26 /20b 04 

0.26/3b b 04 

0.26 752b 04 
0.26768b 04 
0.267846 04 
_ JO . 2 6 8 0 0 b_ 4 _ . 


WibSTaM COORUliAfbS ANO VbLOC I f Y COMPONENTS 
-7.6 7 721b 07 -0.338326 06 0. 28609b 0 3 -Q.754S0E 0j4 

-0 ", 6/7166 07 -0.35032b 06 '0.29961b 03 -0.75459E 04 

— 0 . 67 /llb 0 / -0.5623 2b 00 ■ 0.31313b 03 -0.73438b 04 

-m . o 7706t 0/ -0.37431b 06 0.32665b 03 -0.75417E 04 

-2.6/ 703, t J)7 -0.38630b _0 6 0. 34 016b 03 -0.753956 04 

—/•. 0/6961: 07 -0,49829b 06 P. 33367b 03 -0.75373b 04 

-4 .6 7690b w? — 0 . 4102 7b 06 0.3a719t 03 - 0 . 7535 IE _0_4_ 

-0.6,684b Hi -0 . 42225b 06 0.38070b 03 -0.75329E 04 

- 0 ,67678b 07 -0.43422b 06 O'.39420b 03 -0.7 3306E 04 

-1, 6/6/lb 07 -0 , 4 4 b20b 06 0.43771b 03 -0.75283E 04 

-0 . 6 /6o5-_07_-0 . 4b8 16t 0 6 0.4 2 1.21 1_ O.3-.10 lZ 52.63 b_.M. 


RE NDE ZVOUS 


FIRSV, SECOND HORN. K l 6 5 I UN 0. 358162b. 02 0.158996b 02 

RtlATIVb POSIIION VbCmR -U. 250600c. 04 ? ; 89793.yi._04, 

RhLATjVb VtLOC'ITY CUM-'PnEnTS -0.319665b 02 0. 132 915. t 03 ( j 

_ 4 i ‘ i I 9. 
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Section 8 
CONCLUSIONS 

Planning an optimal rendezvous and steering the intercepting vehicle 
is a complex problem of calculus of variations. Even a modern computer takes 
too long for solving such problems on a real time basis. 

For this purpose simplifying assumptions have to be introduced, and 
a trade-off between simplicity of the guidance equations and accuracy of the 
results has to be made. The simplest way, the impulse approximation, turns 
out to be insufficient in accuracy for realistic cases. 

The present approach is very successful in coming up with rather simple 
guidance equations due to the use of Levi-Civita' s variables. The accuracy is 
such that a good closed -loop terminal guidance scheme could take over after 
the coast phase. 

However, when larger and faster on-board computers are available, 
it might be worthwhile to seek more sophisticated guidance schemes which 
allow a more precise and a more economical steering of the interceptor. 
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APPENDIX 

COMPUTER PROGRAMS 


Prior to the program listings there follows a table giving the meaning 
of the most important program variables or the corresponding name in this 
report. The input the simulation program requires can be seen in the flow 
chart of Section 7, and the output is self-explanatory. 


Main Program 


DT 
NS TP 
SCH 
TINT 
UPD 
XO, YO 
XT, YT 
UT, VT 
X(l) 

X(2), X{3) 
X(4), X(5) 


updating time interval 

number of Runge-Kutta steps in one updating time interval 

integration step 

current updating time interval 

TRUE during updating period 

initial guesses for s, a 

Cartesian coordinates of the target 

velocity components of the target 

time 

Cartesian coordinates of the interceptor 
velocity components of the interceptor 


AO 

a o 

DO 

D o 

R0 

r o 

A2 

a 2 

D2 

D 2 

R2 

r 2 

A10 

a 10 

DAI 

Aa j 

R3 

r 3 

A20 

a 20 

DA 2 

4a 2 

TAR 

*”tar 

A13 

a 13 

DEL 

6 

TAU 

T 

A23 

a 23 

DG1 

A7 i 

U12 

u 12 

BIO 

13 10 

DG2 

A7 Z 

U22 

U 22 

B20 

P 20 

E 

e 

VO 

v 0 

B13 

Pl3 

GAO 

y o 

V3 

v 3 


A - 1 
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B23 

13 23 

GA3 

r 3 

W12 

w 12 

CHO 

*0 

MU 


W22 

W 22 

CH2 

*2 

OMO 

' % 

Z 

Z 

COA 

t 

OM2 

£0— 




coa 


2 



D 

D 

PHI 

9 



Subprogram CONTRL 





H 

initial step size 




NI 

counts iteration cycles 




TOL 

tolerance for 

stopping the iteration 



YOLD 

old approximation for 

a in secant method 



El 

£ 1 

XO 

s o 

Y0 

*0* 

E2 

h 

XI 

S 1 

Y1 

. J* 



X 2 

S 2 

Y2 

«.U 

a 2* 

Subprogram KEPLER 





A 

a 

DT 

At 

U2 

u 2 

A 1 

“l 

OM 

to 

V 

x 2 

A2 

a 2 

R 

r . 

W1 

W 1 

Bl 

h 

SIG 

a 

W2 

W 2 

B2 

P 2 

U 

k \ 

X 

X 1 



U1 

U 1 

Y 

x 2 


A -2 
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C MAIM PROGRAM FOR RENDfZVOUS SIMULATION 

RFAl MU 
logical upo 

COMMON/PARAM / MU , T AU ,£i COA , CHO , CH2 
COMMON/OTHFRS/ VO , V3 , GAO , GA3 , DEL , PH I 

COHMON/TPANSF/ XO 1 YO 1 AO , A2 » OMO , 0M2 ,R0 , R2 , R3 , D , DO , 02 * 

1 DG1 iDG?,OAl » D A2 » I >ES « DET , T AR » Z » oSO » SS2 , 

3 AIO, A20 >B10, B20 » A 1 3 « A23 ,313i B23 , ‘J 1 2 , U22 , W 1 2 , v.‘2 2 

COs' MON/RKG / X ( 5 ) , DX < 5 ) 

READ (5«1) MU « T AU » E , DT , R3 , V3 » GA3 , DEL , RO , VO « GAO , PH I 

1 FORMAT! 4E1 6.6 ) 

WRI TF(6i7) MU , T AU »EiDT »R3» V3t G A3 .DEL , RO * VOi GAG . PH I 

7 FORMAT! 6H1 INPUT/, 

1 2 Oh mu TAU cr DT , 4c 13* 7/ , 

? 2C H R3 V3 GA3 DEL , 4E13.7/, 

3 2 OH RO VC GAO PHI » 4E 13*7/ ) 

C 

NSTp = 10 

TT = O. 

X(l) = 0 . 

X ! 2 ) = RO -* COS (PH I ) 

X ! 3 ) = RO * SINfPHI ) 

X ! 4 ) = VO * SIN! GAO— Ph I ) 

X ( 5 ) = VC * COS!GAC-PhI> 

AO = 1./(?./P0 - VO*\/0/MU> 

XO = 3.141=926536 * SORT ( AO/MU ) 

A2 = l./(2./R3 - V3* V3/MU ) 
yO = 3.1415926536 * SORT!A2/MU) 

WRITFI6, 8) AOi A2 

8 FORMAT ( ///39H SEtMl-MAjOR AXES Of INTERCEPTOR* TARGET. 2E13.7) 

TINT = DT 

UPD -= .TRUE. 

XT = R3 * COS! DEL) 

YT = R3 * SIN! DEL) 

UT = V3 * SIN! GA3— DEL ) 

VT = V3 * COS ! GA3— DEL ) 

WRITEC6. 3) TT.XT.YT.uT. VT 

3 FORMAT <//69H TARGET TIME. CARTESIAN COORDINATES AND VELO 

1CITY COMPONENTS/. 

2 F13.7. 2 !4X. 2Fi 3.7 ) ) 

60 CALL CONTRL 

IF ITINT .LE. DO) GO TO 61 
TINT = DO 
UPD =. FALSE. 

61 CALL KEPLER! TINT, XT .YT.UT.VT) 

TT = TT + TINT 

WRIT£<6. 3) TT.XT.YT.uT.VT 
WRITE! 6. 10 ) 

10 FORMAT !//69H I NT ERC TIME, CARTESIAN COORDINATES AND VELO 

1CITY COMPONENTS/) 
call INTGRTIO.O) 

WRITE (6, 2) X!1),X!2),X!3).X(4),X(5) 

2 FORMAT! E 13.7, 2 ! 4X , 2E 1 3. 7 ) ) 

SCH = TINT / FLOAT!NSTP) 

DO 62 N= 1 » NSTP 
CALL INTGRT(SCH) 

62 WRITE! 6,2) X ( 1 ) , X ! 2 ) , X ( 3 ) » X ! 4 ) , X ( 5 ) 

RO = SORT ( X ! 2 ) **2 + X ( 3 ) **2 ) 

PHI= AT AN2 ( X ( 3 ) , X ! 2 ) ) 

VO = SORT ! X ! 4 ) + X(5)**2> 


A-3 
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GAO - PH I + 1.^707963068 - AT AN8 ( X ( 5 ) « X ( i ) ) 
R3 = SORT ( XT *XT + YT#yT) 

DEL = ATAN2(YT,XT) 

V3 = SQRT(UT*UT + VT*vT) 

GA3= DEL + 1. 5707963268 - AT AN2 (VTtUT) 

IF (UPD) GO TO 60 
WRI TF( 6« 11 > 

1 1 FORMAT < /18H END OF FIRST BURN 1 

CALL KEPLER (C0A+02t XT, YT,UT,VT) 

TT = TT + COA + 02 
„ WR I T E ( 6 i 3) TT,XT,YT,ijT,VT 
WRITr(6, 1?) 


1? FORMAT ( /l tH PFNDF7VOUS) 

CALL KEPLER C COA ,X<2),X(3),XC4),XC5)) 

do = x< n 

x< 1 ) = X ( l > + COA 
CHO = CHp 
WR I TE ( 6' 13) 

13 FORMAT ( /// 1 2H SECOND BURN) 

WR I TE ( 6, 10) ' 

CALL INTGRT(G.O) 

WRITE (6*2) X( 1 ) , X ( 2 ) *X(3)fX(4)»X(5) 


63 


1 4 


SCH = D2/FL0AT (NSTP) 

DO 63 N =1 .NSTP 
CALL INTGRT(SCH) 

WR I TE C 6 , 2 ) XC 1 )«X(2) <X(3)*X(4)iX(5) 

WR I TE ( 6 » 1 2 ) 

TE2 = X < 2 ) —XT 
TE3 = X ( 3 ) — YT 
TF4 = X ( 4 ) —UT 
TEE = X < 5 ) — VT 

WR T TE ( 6 * 14) DO , D2 »X( 1 ) , TE2 » TE3, T E4 « TES 
FORMAT ( /V//29H FIRST, SECOND BURN, MISSION , 
L 29H RELATIVE POSITION VECTOR , 

2 29H RELATIVE VELOCITY COMPONENTS, 

STOP 


3E13 
2E13 
2E1 3 


FND 


7/, 

7/, 

7) 


A-4 
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SUf 'ROUT I NC COM TRL 
REAL VU 
LOGICAL FLAG 

COM^OM/PAPAM / MU*TALJ,E.COA.CHCiCH2 
COM V. ON/OTH F RS / VO. V3»GAG.GA3*OEL*PHI 

COSXON/TRANSF/ XC .YO, AO , A 2 * OMO * CM 2 »ROi i‘- ^ . K3 * L> « UO . u2 » 

DG 1 * DGg iDAl, 0A2 * DLL * DLT , 1 AR * Z * cC>0 .*,.:£* 
AlO*A2O*l3lU,i;.20* A13*A23.*B 13*823, GIZMO.?.: , ■ 1Z, 22 


A >. 

=T 

1 */ ( 2./RC' - VOX- 

O'-'O 

55 

SORT ( .28-XMU/AO ) 

A 2 

= 

1 . / ( 2./R3 - V3* 

'•‘i 

v 


SCRT ( •25*. V «U/ A2 ) 

DE2 

= 

. 3 -X PH I 

TEMP 

= 

SORT ( PC ) 

A 1 0 

“ 

TEMP X COS (DF2) 

A 20 

“ 

TEMP » S1MCDF2) 

OF 2 

= 

GAO - OF 2 

Tr : :p 

- 

. fxtempxvo/O'oc 

C! I W 

= 

i CMP -X SINCDE2) 

P20 

= 

TE-'-’P x- COS (0E2 > 

OFF 

_ 

DEL/2. 

TL-'.P 


SORT ( R3 ) 

A 1 3 

- 

TE'/P * COS (DC 2) 

A 23 

= 

‘TEMP X- SIN(0F2) 

OE2 

= 

PA o - OFF 

tem^ 


,'S X- T p MP X V 3 / CM 2 

'M3 

r 

TF(.’D -X S IN U)F 2) 

P?3 

=: 

TEMP x- COS (DP 2) 


H = XC 

fr i - . --Mr 
F2 - .1 e 'L 
TOL= 1 . LCF — G 
WR I T r ( 6 * A) 

FORK AT (X///8H CONTROL/* &X * 1H*i« i 1A3HSIG* i OX * iH/WX* lHO«SX,2nGC • 
t >X « 3HC OA * 7X » d \ ) 02 * 7X * 3H Ci- 3 * Lj X * 3HCH2 ) 

N I = C 
«=TLF>= w 
T-i_AG= .FALL 1 -:. 

Y - Yw 
CALL FCRX'LY) 

TCI = R*j»DG 2 - UlO-XfOET 
TC^ = r>i'»pc-l + P? '‘-xpry 
CHU = AT AM 2 (TCI* TC? ) 

TCI = npWA? - "’l ? !(CF C , 

; C? = R 2*0 A 1 W 2 2 T-DF S 

CH2 - ATAM2(TC1 *TC2) AT AM2 ( 02 2 * U 1 2 I 

>;.‘h I T E ( 6 * CM XL i Y < Z * D * [)(' t CGA * D2 * C HG * CH2 
r ORA 1 AT ( PF 13. r t FF 1 0.4 , 2F8.4 ) 

I = W I 3 

IF(NI .LF. 7) GO TO /^3 
WITK ( 6* 4 A ) 

FORMAT ( 18; i 7 ITERATION STEPS//) 

C-0 TO F'-> 

ZC = Ei ” A^Sf Z) 

IF ( ZC «LT. TOL ) GO TO 30 

IF (FLAG) STEP = F*Z 
DO 7 1 1=1*7 

ZOLD= z 
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yold= y 

Y = YOLO + STEP 
CALL FCT(XC,Y) 

IF (FLAG) GO TO 83 
F = H / (ZOLD-Z) 

FLAG= .true. 

83 IF (Z .EQ. ZOLD) GO To SO 

IF (ABS(Z) .GT. ZC) Go TO 72 

80 Y'O = Y 
CO = D 
GO TO 73 

72 STEP= -Z / ( ZOLD-Z ) * < YOLD-Y > 

71 CONTINUE 

73 XI = XU + h 
- Y = YO + H 

CALL FCKXl.Y) 

STEP= F * Z 
DO 74 1=1 i 7 
ZOLO= Z 
YOLD= Y 

Y = YULU + 8TEP 
LALL HCI 1X1 i Y ) 

lh IZ ttUi ZULU J GU IO Ol 
It- l ABiU) .bit ZC > GO IO /b 

81 Y 1 = Y 

Cl = D 
GU I U fb 

1 3 8IEH = -Z ✓ IZOLU-ZJift YULU-Y) 

74 CONTINUE 

76 X2 = XG - H 

Y = 2 • *Y’.' - Y 1 

C-^LL HC! ( X2 * Y ) 

S I EP= K * Z 
DU 7 / 1=1*/ 

ZULl'= Z 

YOLD= Y 

■ Y = YOLD + STEP 
CALL FCT(X2.Y) 

I h (Z . EQ. ZOLD) GO 10 82 
IF IARSIZ) .Gl. ZC) GO 10 78 
'82 Y2 = Y 

C? '= D 

GO TO 79 

78 STEP= -Z / ( ZOLD-Z >*( YOLD- Y) 

n con r inue 

79 bIH = -.5*(L1-C2) / ( C 1-2.*C0+C2) 

XU = XU + STP*H 

YU = YU + .b*blP-M Y Jt-Y2+3TP*( Y 1 -2 . *Y0+Y2 ) ) 
H = E2 *STP*H 
GO TO 70 
bO YU = Y 
OE) URN 

end 
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SUBROUTINE FCTCS.SIG) 

REAL MU » LAM 

COMMON/PARAM / MU « T AU , E ♦ COA » CHO , CH2 

COMMON/TRANSF/ XO,YO, AO, A2»OMO,OM2*RO.R2,R3»D,DO »D2* 

1 DG1 ,DG2‘0A1 »DA2,l>ES,DET»TAR,Z,SS0tSS2» 

2 A 1 0 , A20 * B 1 O , 820 » A 1 3 , A23 * B 1 3 * B23« U 1 2 ♦ U22 , W 1 2 « W22 
CO = COS ( OMO*S ) 

SO ss SIN(OMO*S) 

C2 = C0S(0M2*SIG) 

S2 = SINK 0M2*SIG) 

U12 ='C2*A13 + S2*B13 
U22 = C2*A23 4- S2*B23 
W.l 2 =-S2*A13 + C2*B13 
W22 =-S2*A23 + C2*B23 
A 1 2 = CO-s-U 1 2 - S0*W12 
A22 = C0*U22 - S0*W22 
.B12 = S0-K-U12 + CO*W 1 2 
B22 = S0*U22 + C0*W22 
R2 = U 1 2*U 1 2 + U22*U?2 
DAI = A 1 2 - AlO 
DA2 = A22 - A20 
DGI » DA1*C0 + (B12-Bi0)*S0 
DG2 = DA2*C0 + CB22-B20)#SO 
DFT = DGl*B20 - DG2*BlO 
DES = DA 1*W2.2 - DA2*W12 

GO = OM0**3 * SQRT<R0*<DG1*DG1+DG2*DG2) + ( 1 . +AO/RO ) *DET*DET ) 

G2 = 0M2*"*3 * SQRT(R2*<DA1*DA1+DA2*DA2) + U . +A2/R2 ) *DES*DES ) 

LAM = . 25-*MU-*E*S0 
TEMP= LAM + GO/2. 

OO = TAU*GO / TEMP 

D = TAU * (LAM/(LAM+.5*G2) > * C < G0+G2 > /TEMP ) 

D2 = D - DO ; 

B I 1 = B 1 2 + DAl*CO/SO 

B21 = R22 + DA2*C0/S0 
A 1 = (R0+B1 l*BU+B21*B2l)/2. 

OM1 = SORT ( .25#MU/A 1 ) 

SSO = DO/RO 

SS2 = ( 0M2*D2)/<0M0*R2) 

TEMP= 2.*0M1*( S-SS2 ) 

TEN 1 = 2»*OMl*SSO 

COA = OMO/OM1 * (Al*( S-SS0-SS2) 

1 + ( < A 1 0*0 1 1 +A20-SB2 1 ) -K- < COS < TEM 1 )— COS ( TEMP ) > 

2 + { RO— A 1 ) * <SIN( TEMP) -S INK TEM1 ) ) ) / (2.*OMD) 

TAR = A2*S I G + ( < A1 3*B13+A23*B23 )*S2 + (R3-A2)*C2> * S2/OM2 

2 = COA + D - TAR 

RETURN 

FND 
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SUBROUTINE INTGRT (H) 

crkg this subroutine integrates from points y<t> to the point y<t+h>. 

C. THE DERIVATIVES ARE ALSO COMPUTED AT T+H . THE ROUTINE REQUIRES 

C THAT Ni Y. DY BE IN COMMON . AND THAT Y AND DY BE DIMENSIONED 

C BY N. TO INITIALIZE CALL WITH H = 0. 

C ENTERING THIS ROUTINE WITH H=0 CAUSES THE CONSTANTS TO BE SET* 

C THE J VECTOR SET TO ZERO AND THE DERIVATIVES RECOMPUTED VIA THE 

C SUBROUTINE CALLED DEQ 

C COMMON REGION 

COMMON/RKG / X (5) i DX( 5 ) 

DIMENSION A ( 4 ) . BC4). C(4). Q(33) 

IF ( H .EQ. 0.0) GO TO 4 

1 DO 3 I = 1.4 
DO 2 J= 1 . 5 

RRR= A ( I )*DX( J)-B( I )*Q( J) 

Q( J)=0( J) +3.0*RRR-C C I ) *DX ( J ) 

2 X( J) =X< J)+H-*RRR 

3 CALL DEO 
RETURN 

4 D= ■ 7071 0678 1 E+00 
A ( 1 ) = 0.5 

A ( 2 ) = 1.0-D 
A ( 3 ) = 1 . 0+D 

A ( 4 ) = 1.0 / 6.0 

B ( 1 ) = 1.0 

8(2)= A( 2) 

8(3)= A ( 3 ) 

B ( 4 ) = 1.0 / 3.0 

C( 1 ) = AC 1 ) 

C ( 2 ) = A ( 2 ) 

C ( 3 ) = .A<3>. 

C<4)= C.5 
DO 5 J= 1 . 32 

5 0( J)= 0.0 
CALL DEO 
DX ( 1) = 1 . 0 
RETURN 
END- 
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SUBROUT INF DEQ 
PEAL mu 

COMMON/PARAM / MU . T AU ,E» COA » CHO * CH2 
COMMON /RKG / X ( 5 ) ♦ DX ( 5 > 

R = SORT < X ( 2 ) **Z + X(3)**2> 

R = -MU/(P*R*R) 

S = TAU - X( 1 ) 

IF C X C 1 ) .GT. 500.) s = S-C.OA 
S = F / S 
OX ( 2 ) = X ( 4 ) 

DX { 3 ) - X ( 5 ) 

DX ( 4 ) = P*X<2> + S-*COs(ChO> 

DXC5) = P*X(3) + S-^SlMfCHO) 

RETURN 

FND 


SUBROUTINE KEPLER (DT .X.Y.U. V) 

X.Y.U.V ARE THE CARTESIAN COORDINATES AND VELOCITY COMPONENTS OF 
THE VEHICLE. THEY ARE INPUT AS WELL AS OUTPUT PARAMETERS. DT IS 
THE TIME INCREMENT DURING THE KEPLERIAN MOTION. THE GRAVITATIONAL 
PARAMETER MU OF THE CENTRAL BODY MUST BE IN COMMON. 

REAL MU 

COMMON/PARAM / MU . T AU , E * COA ♦ CHO 
R = SORT { X*X+YXY ) 

A = 1. / (2./R - <U*U+V*V >/MU) 

CM = SORT ( . 25*MU/A) 

A 1 = SORT ( • 5* ( R+ABS < X ) ) ) 

A2 = .5*Y/A1 

IF (X .GT . .0 ) GO TO 1 1 

B = A1 

A 1 = A2 

A2 = B 

B 1 = .5*CAl*U + A2*V)/0M 

B2 = .5*(A1*V - A2*U)/0M 

QC = R-A 

QS = A1*B1 + A2-X-B2 

SIG= (DT - .5*QS/OM)/a 
DO 12 1=1.7 

S = S I N C OM*S I G ) 

C = COS ( OM*S I G ) 

F = A«-SIG + S* ( QC*C+QS*S ) /OM - DT 
FP = A + GC*( C*C-S*-S) + 2.-*QS*C*S 
SG = F / FP 

IF (ABS(SG) «LT. 1.5E-I0) GO TO 13 

SIG= SIG - SG 

U1 = A1*C + B1*S 

U2 = A2*C + B2*S 

W1 = R1*C - A1*S 

W2 = R2*C - A2*S 

X = U1*U1 - U2*U2 

Y = 2.*U1*U2 

B = OM / (,5*X+U2*U2) 

U = B*<U1*W1 - U2*W2 ) 

V = R*CU2*W1 + U1-*W2) 

RETURN 

END 


A-9 



